home *** CD-ROM | disk | FTP | other *** search
/ IRIX 6.2 Development Libraries / SGI IRIX 6.2 Development Libraries.iso / dist / complib.idb / usr / share / catman / p_man / cat3 / complib / dlapll.z / dlapll
Text File  |  1996-03-14  |  2KB  |  67 lines
  1.  
  2.  
  3.  
  4. DDDDLLLLAAAAPPPPLLLLLLLL((((3333FFFF))))                                                          DDDDLLLLAAAAPPPPLLLLLLLL((((3333FFFF))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      DLAPLL - two column vectors X and Y, let   A = ( X Y )
  10.  
  11. SSSSYYYYNNNNOOOOPPPPSSSSIIIISSSS
  12.      SUBROUTINE DLAPLL( N, X, INCX, Y, INCY, SSMIN )
  13.  
  14.          INTEGER        INCX, INCY, N
  15.  
  16.          DOUBLE         PRECISION SSMIN
  17.  
  18.          DOUBLE         PRECISION X( * ), Y( * )
  19.  
  20. PPPPUUUURRRRPPPPOOOOSSSSEEEE
  21.      Given two column vectors X and Y, let
  22.  
  23.      The subroutine first computes the QR factorization of A = Q*R, and then
  24.      computes the SVD of the 2-by-2 upper triangular matrix R.  The smaller
  25.      singular value of R is returned in SSMIN, which is used as the
  26.      measurement of the linear dependency of the vectors X and Y.
  27.  
  28.  
  29. AAAARRRRGGGGUUUUMMMMEEEENNNNTTTTSSSS
  30.      N       (input) INTEGER
  31.              The length of the vectors X and Y.
  32.  
  33.      X       (input/output) DOUBLE PRECISION array,
  34.              dimension (1+(N-1)*INCX) On entry, X contains the N-vector X.  On
  35.              exit, X is overwritten.
  36.  
  37.      INCX    (input) INTEGER
  38.              The increment between successive elements of X. INCX > 0.
  39.  
  40.      Y       (input/output) DOUBLE PRECISION array,
  41.              dimension (1+(N-1)*INCY) On entry, Y contains the N-vector Y.  On
  42.              exit, Y is overwritten.
  43.  
  44.      INCY    (input) INTEGER
  45.              The increment between successive elements of Y. INCY > 0.
  46.  
  47.      SSMIN   (output) DOUBLE PRECISION
  48.              The smallest singular value of the N-by-2 matrix A = ( X Y ).
  49.  
  50.  
  51.  
  52.  
  53.  
  54.  
  55.  
  56.  
  57.  
  58.  
  59.  
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.